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Abstract. 

Novae are cataclysmic variable binary systems in which a white dwarf (WD) pri- 
mary is accreting material from a low mass companion. The importance of this accre- 
tion takes on added significance if the WD can increase its mass to reach the Chan- 
drasekhar limit thus exploding as a Type la supernova. In this study we accrete material 
of Solar composition onto carbon/oxygen (CO) WDs of 0.70, 1.00 and 1.35Mq with 
accretion rates from 1 .6 x 10~ 1() to 1 .6 x 10~ 6 Mq yr _1 . We have utilized the MESA stel- 
lar evolution code for our modeling and evolve them for many nova cycles or, in some 
cases, evolution to a red giant stage. Differing behaviors occur as a function of both the 
WD mass and the accretion rate. For the lower WD masses, the models undergo recur- 
rent hydrogen flashes at low accretion rates; for higher accretion rates, steady-burning 
of hydrogen occurs and eventually gives way to recurrent hydrogen flashes. At the 
highest accretion rates, these models go through a steady-burning phase but eventually 
transition into red giants. For the highest WD mass recurrent hydrogen flashes occur at 
lower accretion rates but for higher rates the models exhibit steady -burning interspersed 
with helium flashes. We find that for all our models that undergo recurrent hydrogen 
flashes, as well as the steady-burning models that exhibit helium flashes, the mass of the 
WD continues to grow toward the Chandrasekhar limit. These results suggest that the 
accretion of Solar abundance material onto CO WD's in cataclysmic variable systems, 
the single degenerate scenario, is a viable channel for progenitors of Type la supernova 
explosions. 



1. Introduction 

Cataclysmic variables (CV), including classical and recurrent novae, occur when the 
accretion of material onto a white dwarf (WD), from a companion star, undergoes a 
thermonuclear explosion. It has long been thought that such a process could lead to 
the WD growing to the Chandrasekhar limit resulting in a Type la supernova explosion. 
This is designated the single degenerate scenario. The competing model is the double 
degenerate scenario that involves the merger or collision of two WDs. 

In this paper we investigate the consequence of a WD accreting Solar material 
from a secondary donor star. Many types of close binary systems with a WD primary 
have been suggested as the progenitors for Type la supernovae; thus we have modeled 
accretion onto a wide range of WD masses with a wide range of accretion rates. We 
followed the long-term evolution utilizing the MESA stellar evolution code. In §2 the 
grid of models is discussed and §3 describes the results. In §4 we contrast our results 
with previous work and finish with our conclusions. 
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2. Modeling 



Our mode ling was performed w ith the Modules for Experiments in Stellar Evolution 
(MESA. TPaxton et al]|201lL 12013 ) computer code (version 4798). In our study we used 
initial WD masses of 0.70, 1.00 and 1.35M©. All these WDs were initially bare CO 
cores (C - 0.357, O - 0.619) prior to accretion commencing. The stalling models were at 
a Solar luminosity. Our mass accretion rates were chosen to be 1.6 x 10~ 10 , 1.6 x 10 -9 , 
1.6 x 10 -8 and 1.6 x 10~ 7 M© yr -1 . Additionally, other accretion rates were used in 
order to separate differ ent regime be havior when needed. The material being accreted 
was a Solar mixture dLoddersll2003l Mass loss was triggered when a model reached 
super-Eddington lumin osity. The excess luminosity over Eddington determined the 
rate of mass loss (see IShavivll2002l: fPenissen kov et al.ll2013 ). All our models were run 
for many nova cycles or the time required for definitive long-term behavior to become 
evident. 



3. Results 

In Figure[T](top panel) we show a plot of the surface luminosity versus time for a model 
exhibiting recurrent hydrogen flashes. In the bottom panel we show the mass of the WD 
over the same time period. After the initial growth to the first thermonuclear runaway, 
the model quickly settles into a recurring pattern of flashes in between which mass 
is accreted. During a flash some of t he mass is lost via a super-Eddington wind (see 
Shaviv 120021 : 1 Denissen kov et al.ll2013l) . The positive slope shows that the WD grows in 



mass at a rate of 2.96 x 10~ 8 Mq yr -1 for an accretion rate of 1.6 x 10 _7 Mq yr' 1 . This 
represents an efficiency (mass growth versus mass accreted - mass ejected) per cycle of 
approximately twenty percent. 

In Figure |2] we show a model/accretion rate combination that exhibits an initial 
flash and then settles into an extended period of "steady-burning" where the luminosity 
is constant. This steady-burning continues until either recurrent flashes, as in Figure Q] 
occur, or for higher accretion rates, the model becomes a red giant at which time the 
evolution is halted. During steady-burning the mass of the WD continues to grow as is 
shown in the bottom panel of the figure. 

The 1.35M© WD models exhibit a different behavior at the highest accretion rates. 
After an initial hydrogen flash the model transitions into a steady-burning phase inter- 
rupted by regular helium flashes. During these helium flashes, about half the accreted 
mass is ejected from the WD. However, the WD continues to grow in mass. 

In Figure [3] we plot the growth rate versus the accretion rate for all our models 
that did not terminate in a red giant. The growth rates are much larger for the higher 
accretion rates at a given WD mass. For a given accretion rate, the smaller WD masses 
grow in mass at a greater rate and thus at a higher efficiency than the more massive WD 
models. The 0.70 and I.OOMq WDs, at accretion rates greater than approx. 1 x 10 _7 Mq 
yr" 1 , eventually become red giants at which point we end the calculation. However, the 
1.35Mq models, although exhibiting helium flashes, continue their mass growth at the 
higher accretion rates and never enter a red giant phase. 

In Figure |4] we plot the parameter space we investigated in terms of WD mass and 
accretion rate. We identify the different regimes of behavior in our models and show 
that, for a large range of accretion rates, our models grow in mass and that it is possible 
for a model to start off as low as O.7OM0 and, given sufficient time, will reach 1.35M0 
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Figure 1. Luminosity evolution and mass growth for a 1.35M© model accreting 
at 1.6 x 10~ 7 Mq yr~' exhibiting recurrent hydrogen flashes. The evident decrease 
in mass during each fl ash is caused by mass lost as the outermost layers exceed the 
Eddington luminosity dShavivll2002l) 



and beyond. For higher accretion rates, our 0.70 and 1. 00 Mq models enter a red giant 
phase and cease accreting. The 1.35Mq models, however, never enter such a phase and 
despite recurrent helium flashes their growth continues unabated. 



4. Discussion 



The different regimes we identify in the accretion of Solar material onto a WD are 
summarized in Figure @] The plot has three main regions. At low accretion rates, the 
models, for all masses, undergo recurrent hydrogen flashes. At the highest accretion 
rates for 0.70 and I.OOMq, the models, after an initial flash enter a regime of steady- 
burning which eventually leads to the model becoming a red giant. In such a case the 
accretion ceases and we end the evolution. The third regime, intermediate between the 
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Figure 2. Luminosity evolution and mass growth for a 0.70Mq model accreting 
at 1.6 x 1(T 7 Mq yr" 1 . There are a few hydrogen flashes followed by steady-burning. 



prior two, is where we see steady43urning but the model eventually enters a regime of 
recurrent hydrogen flashes characteristic of a somewhat lower accretion rate where the 
flashes occur from the onset of accretion. 

The intermediate behavior is approximately bounded by accretion r ates of 5 x 10 



to 1 x 10 1 Mq yr . This is similar to, altho ugh slight lower than, th e iNomoto et al. 



d2007h results and the theoretical predictions of lShen & Bildstenl (120071) . For the 1.35M© 
models, we see a similar regime of steady burning that is interrupted by regular helium 
flashes. This occurs only for accretion rates above 5 x 1O" 7 M0 yr" 1 and up to the 
maximum accretion rate we studied of 3. 2 x 10 Mp yr" 1 . 

We find broadly similar results with ICassisi et al. ( 19981) but note some important 
differences. In their plot highlighting the different regimes (Figure 10) they identified 
three main regimes. They found, above a certain accretion rate, that a red giant phase 
ensued without an earlier period of steady43urning as in our results. At low accretion 
rates, they observed recurrent hydrogen flashes that were more powerful at the lower 
accretion rate, in agreement with our work. The intermediate steady43urning regime 
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Figure 3. Growth rate versus the accretion rate for 0.70Mo(solid), 
1 .OOM (dotted) and 1 .35M (dashed). 



they found, however, shrank in size as a function of WD mass until by I.OOMq it 
ceased to exist. We find that this region is approximately constant in size and in fact 
increases in the higher mass models, with helium flashes, so that the red giant regime 
is absent up to the upper limit of accretion rate in our study. I t should be noted that 
they only investigated lower mass WDs and included data from iLivio et alj (119891) for 

I.OO Mq WDs. 

ICassisi et alj (1 19981) concluded that, for the high accretion rate systems, the red 
giant phase resulted in a common envelope that would eject most if not all of the matter 
donated by the secondary star. Thus these systems were not progenitors of Type la 
supernovae. They also concluded that at the lowest accretion rates, responsible for the 
strongest hydrogen flashes, produced a classical novae explosion and that these systems 
could not be progenitors of Type la supernovae. For the models between these limits 
they were limited (by their computing capability) to follow only a few hydrogen flashes. 
It should be noted, however, that for one particular case they followed a model until a 
large helium flash occurred. They surmised that such a helium flash would curtail the 
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Figure 4. Accretion rate - WD mass plane for the range of models we investi- 
gated. The symbols indicate models that become red giants (squares), steady-burning 
followed by hydrogen flashes (stars), steady-burning interrupted by helium flashes 
(triangles) and recurrent hydrogen flashes (diamonds). Between the dotted lines is 
the steady-burning region which does not exist at 1.35Mq. 



growth of the WD and thus such systems were also precluded as being the progenitors 
of Type la supernovae. 

We found no helium flashes in our low mass models for any accretion rate. For the 
steady-burning models they found helium flashes occurred both for low mass models 
(0.516 and 0.6Mq) and the higher mass O.SMq model. They stated that such flashes 
would result in envelope expansion beyond the Roche lobe of the WD and that the 
envelope is lost from the system in a similar manner to the WD becoming a red giant. 

In contrast, the only helium flashes we find are for steady-burning models at 
1.35Mq- One particular case was for a 1.35M0 WD accreting at a rate of 6.4 x 10" 7 Mq 
yr _1 . After approximately seventy-five years of steady -burning the system underwent 
a helium flash in which 2.5 x 10~ 5 Mq was ejected. The radius of the WD peaked at 
0.056/?© an d the peak surface velocity was 10 km s _1 . If this WD were in a CV sys- 
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tem with a typical companion that had a mass of 0.7 'M@ and a radius of 0J5Rq the 
secondary would fill its Roche lobe if the system had a semi-major axis of 2.32/?© an d 
period of 6.9 hours. This gives a Roche lobe radius for the WD primary of 1.01/?©- 
However, our WD model never exceeded more than about six percent of this during 
the helium flash. We also note that the peak velocity of 10 km s" 1 is only about l/30th 
of the escape velocity of the WD which is 2700 km s -1 . We conclude such helium 
flashes do not curtail the growth of our models due to Roche lobe overflow or complete 
envelope ejection. 

The reason for our WD helium flashes being less powerful is that we model our 
mass loss as a super-Eddington wind which prevents the atmosphere from becoming 
extended. Thus a helium flash is not th e dynamic m ass ejection that occurs using a 
mass loss prescription like IPrialnik et al.l (I1978L 1 1979b . In such simulations very large 
radii result from an ini tial thermonuclear runaway and further evolution is curtailed (see 



Starrfield et alj|2012al|b)). More sophist icated treatments of optically thick wind mass 



loss (see Kovetzlll998l : lldan et al.ll2012l ) tend to produce greater mass loss for a given 
accretion rate an d thus correspondingly less efficient growth of the WD primary. In 
Idan et all d2012 ) modeling, steady-burning is quenched and is replaced by a series of 



recurrent hydrogen flashes, which build up a helium layer until a large helium flash oc- 
curs removing all of the accreted mass. However, at lower accretion rates their models 
grow in mass just like ours though less efficiently. 

From our Figure [3]we can estimate the timescale for the growth of systems to reach 
the Chandrasekhar mass and thus be candidate progenitors for Type la supernovae. We 
discard the models that become red giants and the lowest accretion rate models with the 
most powerful hydrogen flashes (accretion rate less than 10~ 8 Mq yr -1 ), assuming that 
they mix in some fashion and become classical novae. The time for the O.7OM0 model 
to attain the Chandrasekhar limit is 8 x 10 6 to 4 x 10 8 years. The LOOM© model takes 
6 x 10 6 to 3 x 10 8 years and the 1.35M0 case requires 3 x 10 6 to 1 x 10 8 years. The 
shortest and longest timescales are seen to vary by two orders of magnitude. 



5. Summary 

We conclude that the single degenerate scenario for Type la supernova production is 
viable as long as the accreted material does not mix with the core material. Not all 
combinations of WD mass/accretion rates lead to continued mass growth of the WD 
primary but some do allow accretion to continue unabated. However, at lower WD 
masses, high accretion rates lead to red giant formation of the WD, which through 
Roche lobe overflow and a subsequent common envelope phase eject mass from the 
system. Though we did not investigate very low accretion rates, where diffusion mixing 
becomes important, it is suspected that the more powerful hydrogen flashes would eject 
enough mass to prevent growth to the Chandrasekhar limit. 

Our models that undergo recurrent hydrogen flashes or steady -burning followed 
by either hydrogen or helium flashes continue to grow in mass. We find that the he- 
lium flashes are not po werful enough to curtail the mass growth. The first point is in 
agreement with that of lldan et all ( 2012 ) though the second point, in regard to helium 



flashes, is at odds with their work. This is most likely due to the different approach to 
mass loss; their optically thick wind model as opposed to our super-Eddington mass 
loss prescription. 
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We find that for our WD models that grow in mass; the timescale for reaching the 
Chandrasekhar limit varies by two orders of magnitude, dependent upon initial mass 
and accretion rate, from 3 x 10 to 4 x 10 8 years. These results suggest that Solar type 
material accreted onto a CO WD is a channel for Type la supernovae production. 
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